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Short thermalization times of less than 1 fm/c for quark and gluon matter have been suggested by 
recent experiments at the Relativistic Heavy Ion Collider (RHIC). It has been difficult to justify this 
rapid thermalization in first-principle calculations based on perturbation theory or the color glass 
condensate picture. Here, we address the related question of the decoherence of the gluon field, which 
is a necessary component of thermalization. We present a simplified leading-order computation of 
the decoherence time of a gluon ensemble subject to an incoming fiux of Weizsacker- Williams gluons. 
We also discuss the entropy produced during the decoherence process and its relation to the entropy 
in the final state which has been measured experimentally. 



I. INTRODUCTION 

Collisions of nuclei at very high energies have been 
studied at the Relativistic Heavy Ion Collider (RHIC) 
in recent years, in order to explore the formation and 
the properties of the quark gluon plasma (QGP). One 
striking discovery was the fact that ideal hydrodynam- 
ics could describe many salient features of the expansion 
and cooling of the fireball [ll, S S]- ^^ particular the 
azimuthal asymmetry in collisions with nonzero impact 
parameter requires an onset of the hydrodynamic expan- 
sion at rather early times, between 0.5 and 1 fm/c. This 
represents a puzzle since the application of ideal hydrody- 
namics mandates complete thermalization of the system. 

Scenarios based on the picture of a dilute system of 
perturbatively interacting partons and minijets fail to 
describe early thermalization [4] and the formation of 
anisotropic transverse flow [5,] . Multi-parton interactions, 
either involving ternary collisions [y, Ml or collective ef- 
fects mediated by plasma instabilities [S] have been pro- 
posed as solutions to this problem. However, both mech- 
anisms require a large initial entropy production before 
they can become effective. The color glass condensate 
model [3, llO] seemingly offers a viable explanation of the 
conundrum. It introduces the saturation scale Qs which 



sets the scale for all dynamical processes. For RHIC [11 1 
the magnitude of Q^ is estimated to be approximately 
2 GeV^, and thus thermalization times of order \/Qs do 
not seem impossible. However, the production of entropy 
is a nontrivial problem in any model that is based on the 
assumption of the dominance of classical fields in the ini- 
tial state, like the color glass condensate model pj|, |lj] . 
The situation is further complicated by the fact that 
it is not clear how close quark-gluon matter must be to 
complete thermalization for hydrodynamics to be suc- 
cessful. It is also conceivable that a hydrodynamic evo- 
lution starting at a later time and supplemented by other 



mechanisms of transverse dynamics during the long off- 
equilihrium phase, can give an equally good or even bet- 
ter description of the data. For example, the anisotropic 
collective transverse flow may be generated, in part, by 
interactions of minijets with the bulk medium |14| or 
by anisotropics in the initial gluon field [15|. Hydro- 
dynamic calculations with viscous corrections describing 
small deviations from equilibrium start to become avail- 
able [i^, [13, 13 ^'^^ ™^y soon help to test this possibility 
quantitatively. 

One necessary ingredient for thermalization is the de- 
coherence of the initial gluon field. Coherent fields can 
lead to large anisotropics in pressure and even negative 
pressure, which are symptoms of a state very far from 
thermal equilibrium. Thus the decoherence time Tdoc 
should be even smaller than the equilibration time rth. 
We argue that the fundamental process at work is some- 
what analogous to Coulomb explosion imaging, see [l9[, 
used routinely in molecular physics. If a molecule trans- 
verses a very thin metal foil all bonds are broken, the 
ions decohere and fly apart. From the momentum distri- 
bution of the fragments one can then extract information 
on the original wavefunction. In heavy ion collisions each 
Lorentz contracted ion acts like such a foil for the other. 

That the loss of information due to decoherence can 
generate a rapid increase in entropy in early phases of 
heavy ion collisions was realized already early on, e.g. 
by Elze [20|. The fact that the hydrodynamic evolution 
is known to be very close to the ideal one and thus isen- 
tropic in heavy ion collisions, further stresses the need for 
massive entropy production in very early phases. Indeed, 
two of us have argued in [21], that decoherence can easily 
generate a large fraction of the total produced entropy. 
However, the question of the appropriate time scale of 
entropy production remained open. 

A computation of the decoherence time in leading or- 
der in perturbation theory was recently presented by two 



of us in [23] ■ Here, we want to strengthen this argument 
by presenting a calculation of the decoherence time as 
a function of the gluon 2-point function in the nucleus. 
We then proceed to evaluate our general result within the 
framework of the McLerran-Venugopalan model [9,] . This 
is made possible by new results for the effects of the run- 
ning of the coupling constant on the color glass conden- 
sate [2j|, which solve a hitherto unresolved UV problem. 
Our result agrees with that of Ref. [23| within theoretical 
uncertainties and suggests that, indeed, t^cc '^ Qj^j ^-nd 
that Tfjec is numerically smaller than 1 fm/c at RHIC. 
In Sec. V we can then revisit some of our previous argu- 
ments [2IJ about entropy production in the framework of 
the short decoherence times at RHIC. 



II. THE DECOHERENCE TIME 
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FIG. 1: Lowest order perturbative process that contributes 
to the time evolution of the matrix element D^ ^ for fields A, 

A which are separated from the initial fields A' and A' by a 
rapidity gap. 



We describe the gluons in a nucleus by a density oper- 
ator D. This nucleus is subjected to the incident gluon 
field of a second, large and very fast nucleus scattering off 
it. We compute the time evolution of the density opera- 
tor D under the influence of the perturbation presented 
by the second nucleus. In the following we denote the 
initial unperturbed gluon field of the first nucleus with 
A' , the final gluon field with A and the field of the second 
nucleus with B. The decoherence time of the gluon field 
A is defined as the inverse decay time of the ratio 



TrZ^^ft) 



(1) 



Let us introduce some useful notations. We deal with 
matrix elements D^ ^ — {A\D\A) of the density matrix. 
We can treat the final gluon field as almost on-shell for 
long times after the collision. In practice that means 
that we can decompose it in free modes \A) = \k,X,a) 
characterized by momentum A:, polarization A and color 
a. We can use the same technique for the field B of 
the fast moving nucleus 2 which is Weizsacker- Williams- 
like. On the other hand, the initial-state gluons A' in 
nucleus 1 are in a bound state and generally off-shell. Wc 
do not attempt to describe this field in detail. It turns 
out that the only two ingredients needed are an ansatz 
for the matrix elements D^, ^, of the density matrix of 
the bound fields, and the matrix elements Ha',a of the 
Hamiltonian coupling the fields A, B and A' . 

Wc arc interested in processes in which gluon modes A' 
of the nucleus at rest (which are centered around rapidity 
y = 0) are scattered into modes A with large longitudi- 
nal momenta (K > 1) so that the overlap with the initial 
state is very small. In that case the leading contribution 
in the time evolution comes from second-order perturba- 
tion theory 






nucleus 2 



nucleus 1 



/ 




1 




1 








1 


/♦♦k 


'♦k- 


y(2) 




y(i) = o 




X 3 axes 







FIG. 2: Sketch of the relevant rapidity distributions. The 
rapidity of the gluon 1 in the nucleus at rest fc'(l) is close to 
zero. The Weizsacker- Williams gluon with momentum I has a 
nearly boost-invariant rapidity distribution. The distribution 
of the final gluon momentum fc(l) is, therefore, also nearly 
boost-invariant. 



X H,^,{i)D^,,,{0)HA'At') (2) 



Note that we have suppressed the field B in the notation 
of the matrix elements. We treat the field B rather as 
an external parameter given by the second nucleus. The 
interpretation of this process is illustrated in Fig. [TJ The 
rapidity distributions of the initial and final gluons are 
schematically shown in Fig. [51 



Wc have to specify the relevant part of the Hamiltonian 
matrix element Ha',a which is the three-gluon vertex 



with fields A, A' and B. It is given by 



(2^)4 V2fc+1^ 



x/-X(0 



(fc + fc')^e"*(fc, A)e,(fc', A'K" (fc', A') 
- (fc + iye''*{k, A)e^(fc', A'K'^'lfc', A') 
+ (/ - k'Yetik, A)e^(fc', A')yl'°'''(fc') 



(3) 



Here i" are the adjoint 5'C/(3) generators, /"^'^ are the 
structure constants of SU{3) and 5 is the coupHng con- 
stant. Note that we have used Fourier transformations 
of the operators of the initial gluon field and the field of 
the second nucleus 



(4) 



xej,(fc',A')i"'+h.c.') 
B^i^) -T.I]0jA (e-^'-"B^^(Oi^ + h.c.) (5) 



with operators A'j^, y and B'^^{1). This is similar to the 
usual expansion of free fields, which we apply in our case 
to the final gluon field A. Eq. [3]is then analogous to the 
familiar case of three interacting fields that are asymp- 
totically free. However, in order to respect the unknown 
dynamics of the initial field we do not specify the result 
of the operators A'§ y and B'^^{1) acting on the initial 
states. Rather, we use the matrix elements 



A"'ik',y) = {0\Ai^y\A'), (6) 

S"^(/) = (0|B'='^(/)|B). (7) 



for given initial states \A') and \B) in Eq. [3l 



We now consider measurements made at times much 
larger than the time it takes for the field A' to interact 
with the Lorentz-contracted fast nucleus, i.e. L dt' — > 
dt' . In this case the full integral d'^x over each three 



x 

gluon vertex can be easily carried out. We define the 
evolution matrix for the time evolution of the density 
matrix elements as 



W 



AA'AA' 



dm A' A I dtH\, ^ 



,rt,r 



.2,c.'a,.a'arf^fc' d^k' ^,.',,^, ,,,,ra'\v „^B l{k^k')Bl\k-k')), 

2vVk+k+{VT)^ 



,g2jcaajc 



(27r)4 (27r) 



-^'"(fc',A'K" {k',X')- 



eUk,X)ea{k,X)eUk',X')el{k',X') 



X [{k + k'Yg"" - {2k - k'Yg^t' + {k- 2k'yg^"'] {k + k'fg^^ - (2fc - k'fg'^'^ + (k ~- 2fc')*/' 

I 



(8) 



Since the second nucleus is moving extremely fast and we 
are not interested the time evolution of its gluon fields 
B, we have averaged over the fields in the second nucleus 

B1,{l)Bl'' (i) -^ {Bl{l)Bf{i))^, leaving us with an expres- 
sion that only depends on the initial and final fields in 
nucleus 1. 

Arguing with translational and rotational invariance 
in the transverse plane [24!| and using that it is moving 
along the light cone, we can decompose the two-point 
correlation function of the gluon field of nucleus 2 as 



p^q^ 



5-'^5^^5^j{2t:Y5\ip^~ci^) 

■K'^5{p~)5{q-)p^pj 



p+q+ 



Pi 



G{pi_). (9) 



where G(j)±) is the scalar correlation function for the 
gluon field in the fast moving nucleus, and i,j denote 
the transverse directions. 

Thus far we have not made use of the fact that nucleus 
1 is at rest and the results are valid in general as long as 
the phase space of initial and final gluons is sufficiently 
different. Now we note that the final gluon momenta 
have large -I- components, much larger than the original 
ones, i.e. fc+,fc+ > k'+ ,k'+ with k^ = {k° ± k^)/V2. 
Therefore the dominant terms in Eq. ([5]) are those with 
the maximum number of factors A:"*" or k~^ . These are 

i2k-k'Ye^{k',X') w 2fc+e-(fc',A'), (10) 

(2k~k'fe^4k',X') ~ 2fc+e;(fc',A'). (11) 

Hence, the leading order contribution to the evolution 
operator for the gluon density matrix in our specific kine- 



matic situation is 



2 rca a sca'a " '^ " '^ 



^AA',AA' - .9 / / (27r)4 (27r)4 

xG(lki-kll) '!L_-A"'\k',X')A'^'\k',X') 

2VVk+k+ 

xe-{k',\')e-*Ck',X'), (12) 

where we have introduced the abbreviation 

V .vi. xv^-v^ ^ (kx-kir(kx~kip 

Note that V is essentially the product of the projections 
of the two polarization vectors onto the transverse direc- 
tion given by the vector kj^ — k'^ = k^ — k^ . 

III. THE GLUON DENSITY MATRIX 

The evolution of the gluon density matrix of nucleus 1 
can now be computed through 



DAjt)= E^1'.A'(0)W^AA',AA'W' 



(14) 



A\A' 



For our calculation, we do not need the nuclear den- 
sity matrix itself, but just the expectation value of 

A'1., yA'y y in the ground state of nucleus 1. As in 
[22| we use an ansatz for this expectation value which is 
diagonal and exhibits a Gaussian momentum distribution 
with width l/(^: 



■^ k'.X'^ k'.X' 



Y, D^,^^,{QmA'i,y\A'){A'\Ai:y |0) 



A', A' 



A',A' 



(15) 

In order to determine the normalization constant J\f 
we could calculate the energy density of gluons in the 
nucleus, pg = ^ — y^^ J dtTr[HD], but we will see 
later that our final result does not depend on TV. 

Returning to the expression (fT4|) for the final-state den- 
sity matrix, we obtain 



-Dj,(*)=9'W< 



/I^^Ei'-t^'.^o 



X {2Try S{k- - k'~)S(k- -k'-) 



,_,(52(ki-ki) 



2V^/k+k+ 



x<5a,G(|k^-kl|)AAC'e-';^('="''+''") 



(16) 



The (5-function enforces k' = k . Furthermore, we can 
use k K, k to argue that V w \^x\ because the projec- 
tions in V are maximal if A = A and the average value 
upon integration over the directions of k'^ should be 
(cos^ (p) ~ 1/2. The same argument also allows us to use 
the approximation ^^i |e~(A;', A')p « 3/4, allowing for 
three polarization states of the off-shell gluons in nucleus 
1. We can thus write 



ig^N,C^f 



^^'^ " 327r3/2 



haS^xe 



<\k-f 



.8ik--k-f^'^--^^ 



2V\/k+k+ 



f d^k'^G{\k^-k'j 



(17) 



We now introduce the convolution of the gluon two-point 
function with the Gaussian profile 

F{k^) = Jd^k'^GQk^-k'^De-i'"'^", (18) 

and thus obtain our final expression for the final-state 
gluon density matrix (|14p : 



^A.A ^ ■ 



A^|^^..^a„^^^e-C^(-)V(fc,), 



8^V 



2\/k+k+ 



(19) 



where the three-dimensional delta function refers to the 
"— " and "_L" components of the momenta. We note that 
the density matrix resulting from the interaction with the 
external gluon field for large times is diagonal in all quan- 
tum numbers except for the longitudinal momentum, in 
accordance with the result obtained in ref. [22]. 
Now we proceed to calculate the traces. 



_^.,,..,r 3asN,iN^^l ) f dk+ f d^k^ 
^ ' 32^2 



k+ 



^/(#^(^-)- 



(20) 



The integral over k~^ should be regulated by a phase space 
projection, because the gluons are (almost) on the mass- 
shell, but there is no need to specify the details here. We 
have also assumed that the final-state gluons can carry 
only transverse polarizations, because they are nearly on 
mass-shell. The trace of the square is 

TvD^= Y,Da,A'Da',a 



A,A' 



^T^6^{Q)Af' 



.9alN^iN^-l)C 



128V2^(27r)6 



dk+\^ r d^k± 



{2^y 



F{k^f (21) 



For the ratio that we are seeking this leads to the expres- 
sion. 



TrD'^ _ 1 C\/2^ h 

[TrL>]2 ~ (27r)3(53(o) 2{N^ - 1) C^ 



where 



l2 = j d^k^9^F\ki_) 



(22) 



(23) 
(24) 



We note that the volume associated with 5^(0) = 
(5(0~)52(Ox) is proportional to the light-cone "time" 
variable a;+, which is conjugate to k~ . When boosted 
to mid-rapidity, x"*" transforms as x'+ — x+/(27), un- 
like the regular time coordinate a;*^, which transforms as 
x'^ = ^(if — l3x^). This fact is in line with our intuitive 
picture in analogy to Coulomb-explosion: Decoherence 
occurs due to the fact that Lorentz contracted nucleus 
2, which acts as if consisting of incoherent color-charges, 
passes through nucleus 1. 

The result we obtained coincides with that of our pre- 
vious calculation in predicting a decoherence behavior 
^ l/x~^ at leading order [22]. It is given as a function of 
F(k±) which can be evaluated using different models for 
the initial two-gluon correlator G. 



IV. DECOHERENCE IN THE 
MCLERRAN-VENUGOPALAN MODEL 

In this section we want to compute /i and I2 using the 
standard two-point gluon function from the McLerran- 
Venugopalan model. We work with the Fourier transform 
G(p) = J (Pxe~^^^f(x). Here and in the following we 
suppress the index "_L" for easier notation. / has first 
been calculated in the McLerran-Venugopalan model in 
|25| . We follow the conventions in Lappi ^4] and write 



/(x) 



m'c 



1) 



Ncg^x 






1 



-g'^Nc/{S7r)ti^x'^ In 1/(2;A) 



(25) 



where /x^ is related to the saturation scale Qs ^ g^f^- 
This result is only valid for x < 1/A where A is a IR 
cutoff and we set / = for a; > 1/A. All vectors are 
2-vectors in the transverse plane. 

Note that we have defined factors of the coupling con- 
stant g into /i and I2 without canceling them in the ra- 
tio. We do so because only the square of the gluon field 
strength tensor times the running coupling as has well- 
defined properties and we consider g'^G{p) to be the phys- 
ical quantity. The correct implementation of the running 
coupling is a topic of intense investigations. We follow 
the prescription by Kovchegov and Weigert [2^ and sub- 
stitute 



9 



\K')g\l/x^) 



(26) 



For /i we obtain 



h = (2^)2- lim(gV(x)) 



(2^)^-lim^ 



N' 



X^ K^o 27r 



^i'g\^')g'{l/x')\nl/{xK) 



(27r)V(A^) 



(iV2 - 1)^2 



/3oA2 



where we used the 1-loop running coupling with 



^° = y^-^^- 



(27) 



(28) 



Obviously this is a well defined expression, while /i with- 
out the running coupling would have led to a logarithmic 
UV divergence. 

On the other hand, after two Gaussian integrations we 
see that 



(2^) 



(2^) 



<Px 
X^ J (2^)2 

3^2l6(7V2-l) 



-£C^/(2A^) 4 /•2 



g'n^) 



2 /.A 



A4 



m 



dx 



-x^/{2\^) 



2/- a2\,,2„2 



1 -exp -5^(A^)^^a; 



7riVcln(lAT2A2) 
/3oln(l/:E2A2 



QCD-' 



(29) 



The scale A should be chosen such that for typical trans- 
verse momenta 1/x ideally fulfills [2^ 



1 



> A > Aqcd ■ 



(30) 



For X between and 1/A the ratio of the logarithms 
should hence obey 



0<H1/x'A')/HI/x'AIcd)<1,: 



(31) 



We conclude that I2 < I2 where I2 is given by the last ex- 
pression in Eq. ([29]) with the ratio of logarithms replaced 
by one. This inequality is useful because the remaining 
integral in I2 can be solved analytically. 

After replacing u = x^ !s3' and introducing the short 
notations a = (2A2A2)-i and h = T:N^i?g'^{K^)/{k^l3Q) 
the integral is 



/2 = 2(2^) 



A4 7V2 J^ y2 



X fe-"" - 2e-(''+^)" -f g-(a+26)«\ ^ (32) 
After two partial integrations this can be brought into 



the form 



h = 2(27r) 



,A'{N^-iy 



A4 7V2 
dw In u (a^g-"' - 2(a + 6)2e{'»+6)" 



(33) 



Now we expand the integration region of the remain- 
ing integral to infinity, anticipating that the integrand 
is rapidly vanishing for u —^ oo. Then the integral gives 



2(27r)' 



A4 m 



_ -au I 2g(a+b)i' _ -(a+2b)u 

+ a(7 + lna) - 2(a + 6) (7 + ln(a + 6)) 

+ (a + 26)(7 + ln(a + 26))] (34) 

where 7 is Euler's constant. 

In the region of applicability for the Color Glass Con- 
densate we expect h ^ a ^ C>{1). Hence the bracket 
in the last equation is to good approximation equal to 
2&ln2. Finally we obtain 



/2 w21n2(27r)V(A^) 






(35) 



giving the following bound for the relevant ratio of inte- 
grals: 



< 



2/3oln2 



(/l)2 (2^)2g2(A2)iVe^2 



(36) 



We now return to Eq. (|^ . The expression (27r)^J^(0) 
in the denominator gives the transverse normalization 
area times the observation time T . To fix the transverse 
normalization area one can follow two different lines of 
argument. The incoming giuons from nucleus 1 are ef- 
fectively localized within the transverse area ttC,"^ given 
by the initial density matrix (|15p . If the density matrix 
D is interpreted as that of a completely coherent system 
TrZ?2 « (Tr_D)2 than the transverse normalization area 
has to be chosen as tt^^. If, on the other hand, one prefers 
to extend the normalization area to the whole area of the 
nucleus, 7ri?2, one has to take into account that the start- 
ing value of TrI?2/(TrI?)2 is not close to one but rather 
of the order C^/i?^ and one should thus ask after which 
time i?2xr 1)2/(2 (XrD)2 j^g^g dropped to 1/e. (In the 
latter case D has the form of a block-diagonal matrix 
with B? /C,"^ blocks.) It is reassuring, that both lines of 
argument lead to the same result. 



< 



2/3oln2 



1 



i?2 TvD'^ 

C2"[Tri:i]2 " (27r)5/2iVc(7V2-l),g2(A2)^2(7- 



(37) 



Defining the decoherence time Tdoc as the time where this 
ratio has dropped to a value 1/e and fixing the physical 
saturation scale as Qs = g^{^^)^J,, we obtain the upper 
bound 



Tdc 



< 



8eln2 \ / g^i^^l 

CQs) Qs 




ci) i ■ ^''^ 



All factors in parentheses being of order unity, we thus 
conclude that Tdoc ~ Qj^ in agreement with the result 
obtained in Ref. I2al. 



V. DECOHERENCE ENTROPY 

We now turn to the question how much entropy can 
be produced by the rapid decoherence of the initially 
coherent nuclear gluon field. In order to illustrate the 
mechanism, we first discuss a simple model for which the 
relevant calculations can be performed exactly [2l|, but 
which is sufficiently general to permit a semi-quantitative 
estimate of the entropy produced by decoherence in a 
heavy ion reaction. 

The quantum mechanical analogue of a classical field 
is a coherent state [2a| 



l*[-^]) == nexp(iakAakA ~ i"kAakA)|0), 



(39) 



k,A 



where the amplitude a^x is determined by the classical 
current J creating the field 



QfkA = (^t^kV) ^/^EkA • J(k, UJk.) 



(40) 



Let us begin by considering a single mode kA. The co- 
herent state can be written as a superposition of particle 
number eigenstates: 



a 



aP/2V^ " 



(41) 



Being a pure quantum state, \a) is described by a density 
matrix 



Pmn = {m\a){a\n), 



(42) 



which satisfies the relation p^ — p and has no entropy: 
S = —Trp Inp = 0. 

Complete decoherence of this quantum state corre- 
sponds to the total decay of all off-diagonal matrix ele- 
ments of the density matrix, yielding the diagonal density 
matrix 



dec 

rrnn 



'6r, 



\a\ 



2 a 



2r! 



n! 



(43) 



The particle number in this mixed state follows the Pois- 
son distribution, and the average number of particles is 
fi = |ap. The entropy content of the mixed state is given 
by 



S. 



(cs) 
dec 



n=Q ^ 



°2, n" 
= e"" > — -(nlnn- n- Inn!), (44) 

where the superscript "cs" indicates that the result holds 
for a coherent state. With the help of Stirling's formula 
and the integral representation of the logarithm. 



ds 



(45) 



the sum in (|44|) can be performed yielding an analytical 
result that is valid asymptotically for n 3> 1 (actually, 
the approximation is excellent already for n fn I): 



^(cs) 
^dec 



i(ln(2..)+l-^ 



(46) 



It is not surprising that the entropy is proportional to 
In v^, because we have deleted all information about the 
relative signs of the amplitudes {a\n) by eliminating the 
off-diagonal elements of the density matrix. The num- 
ber of significantly contributing elements is given by the 
width. An = \/E, of the Poisson distribution. That the 
decoherence entropy is controlled by An, rather than by 
n, can be seen by considering more general pure quantum 
states, for which the average occupation number n and 
the occupation number uncertainty An are not related. 
For a pure state with n ^ An 3> 1 and in the Gaussian 
approximation, it is straightforward to show that the de- 
coherence entropy is given by 



Sd 



(ln(2^(An)2) +1 + ...) 



(47) 



confirming our assertion. For a classical coherent state 
(HTjl . the expression (|47p coincides with (P5|) . 

We also note that the entropy for a single quantum 
oscillator in equilibrium at temperature T is given by 



5 = ln(n + 1) + n In 1 + - 
' n 



(48) 



where n = (e'^'^ ~ 1) ^ is the average occupation num- 
ber. Asymptotically, for large n, one obtains ^oq ~ 

2'S'dec ' i- ^- ^^^ thermal entropy becomes twice as large 
as the decoherence entropy. However, for small to mod- 
erate occupation numbers the ratio S^^^ /Scq is close to 
unity. Figure [3] shows the decoherence and equilibrium 
entropies as a function of the average occupation number 
n. For not too large values of n, the decoherence process 
generates a large fraction of the equilibrium entropy, and 
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FIG. 3: Decoherence entropy ^dcc for a coherent state of a 
single field mode and equilibrium entropy Soq for the same 
average total energy as a function of the average occupation 
number n. 



any subsequent equilibration process adds only a small 
amount of entropy to it. 

What does this imply for quantum field theory, where 
the field is a system of infinitely many coupled oscilla- 
tors? Assume that, after decoherence, the system can be 
described as a collection of N particles, given by some dis- 
tribution function over single-particle states, which were 
generated by the decoherence of A^cs coherent quantum 
states. Examples of such states are the internal wave- 
functions of nucleons forming a large nucleus, or a quark 
with its comoving gluon cloud. Each coherent state con- 
tributes on average n = N/Ncs partons. Then, after 
full equilibration, the thermal entropy is of the order of 
Sth ^ Ncsfi = N, while for the decoherence entropy we 
get ^doc ~ ^cs5 ln(27rn). The ratio of the two entropies 
is 



S'dcc ln(27rn) 



5th 



2n 



(49) 



i. e. for large amplitude quantum states, which turn into 
many particles per coherent mode, the decoherence con- 
tribution to the thermal entropy is small. On the other 
hand, if the individual occupation numbers are of order 
one, the contribution is sizable. This case applies to our 
problem of interest, the collision of two nuclei at high 
energy, as we will discuss in the following. 

For the coherent color fields in colliding nuclei, the 
average number of decohering gluons per transverse area 
has been given by [27| 



dN 
d^bdy 



CF\n2Ql 



(50) 



where Cp = 4/3. The characteristic transverse area, over 
which the color fields in nucleus 2 are coherent, is tt/Q^, 



and one can argue that the longitudinal coherence length 
is of the order of Ay « l/a^ @. We thus obtain for 
the average number of decohering partons per coherence 
domain 



dN 



(Pbdy Q2 



Ay; 



Cpln2 



nat 



(51) 



For this value, our arguments presented above indicate 
that the entropy produced in the decoherence process is 
about half of the equilibrium entropy. Applying Eq. (1461) 
and using that the initial number of coherent domains per 
transverse area is (QsR)^, we find that the total entropy 
per unit rapidity produced by decoherence in a Au+Au 
collision at RHIC is 



dSda 

dy 



2z?2 



QiR 
2Ay 
QlR'a 



(ln(27rn) + 1) 
2CFln2 



In- 



1 



at 



1500, (52) 



where we used the values [JtJ Q^ « 2 GeV^, i? = 7 fm, 
and as ~ 0.3. This value accounts for about one-third 
of the entropy measured in the final hadron distribution 

MM- 



VI. CONCLUSIONS 

We advocate the idea that a large fraction of the total 
entropy produced in high-energy heavy-ion collisions is 
generated by decoherence of the many-body quark-gluon 
wave functions of the colliding nuclei in the very first 
phase of the collision. 

We presented an improved determination of the deco- 
herence time Tdec as a function of the initial gluon cor- 
relation function. Within the color-glass-condensate for- 
malism this leads to a decoherence time Tdcc "£ 1 fm/c 
which agrees with the result of an earlier calculation. We 
also estimate the entropy produced through decoherence 
of the initial gluon field and find that it could contribute 
about one third of the total entropy observed at RHIC. 
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